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ABSTRACT 

We present the results of a spectroscopic multisite campaign for the /? Cephei star 
12 (DD) Lacertae. Our study is based on more than thousand high-resolution high 
S/N spectra gathered with 8 different telescopes in a time span of 11 months. In 
addition we make use of numerous archival spectroscopic measurements. We confirm 
10 independent frequencies recently discovered from photometry, as well as harmon- 
ics and combination frequencies. In particular, the SPB-like g-mode with frequency 
0.3428 d~^ reported before is detected in our spectroscopy. We identify the four main 
modes as (Zi,mi) = (1,1), (^2,n^2) = (0,0), (^3,n^3) = (1,0) and {li,m4) = (2,1) 
for /i = 5.178964d-i, /s = 5.334224d-i, /a = 5.066316d-i and fi = 5.490133d-\ 
respectively. Our seismic modelling shows that /2 is likely the radial first overtone and 
that the core overshooting parameter aov is lower than 0.4 local pressure scale heights. 

Key viTords: stars: main-sequence - stars: individual: 12 Lac - stars: oscillations - 
stars: interiors 



1 INTRODUCTION 

12 (DD) Lacertae (12 Lac) is one of the best observed 
(3 Cephei stars, a class of variable early B-type stars 
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whose light and radial velocity changes are due to grav- 
ity and p ressure-mode pulsations of low radial order 
IStankov fc Handler 20051 ). The variability of 12 Lac has 
been known for one century duri ng which the s t ar ha s been 
extensively studied. We refer to [Handler et al.l (|2006h for a 
very detailed summary of the work accomplished in the past. 

Despite the numerous earlier studies, safe mode iden- 
tification of the observed modes was not achieved un- 
til recently. T hanks to an intensiv e photometric multi- 
site campaign, [Handler et al.l (|2006l ) unambiguously iden- 
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Figure 1. The average profiles of the Sim A4553 A line from 
each observatory. The dashed line shows the average profile from 
the BAO observatory. This profile deviates from all the others 
because the relatively small number of spectra implies that the 
beat cycle is not well covered (see Table[TJ. 

tified the ^-values of the five modes with the high- 
est photometric amphtudes. In addition, they found con- 
straints on I for the six other independent modes de- 
tected in their dataset. Their ^'-identifications ruled out 
the assumption previously adopted for stellar modelling 
(e.g. iDziembowski fc Jerzvkiewicj Il999t ) that three of the 
strongest modes, almost equidistant, belong to the same 
multiplet. Indeed, these three modes are actually associated 
to three different values of the degree I. 

Obviously, reliable empirical mode identification was in- 
dispensable before any attempt of in-depth seismic mod- 
elling of 12 Lac. To complement the photometric results, a 
spectroscopic multisite campaign has also been devoted to 
the star. The additional constraints are presented in this pa- 
per. They mainly concern the identification of m-values for 
the strongest modes and the derivation of the stellar equa- 
torial velocity. A det ailed abunda n ce an alysis of 12 Lac was 
already presented in I Morel et al.l (|2006l ) , showing that the 
abundances of all considered chemical elements are indistin- 
guishable from the values found for OB dwarfs in the solar 
neighbourhood. 

Besides our line-profile study, we also describe a detailed 
stellar modelling based on all available observational results, 
state-of-the-art numerical tools and up-to-date physical in- 
puts ap propriate to model /3 Cephei stars, as e xplained in 
detail in lMigliol (|2007l l. IDziembowski fc Pamvatn vkh (2008) 
already computed models for 12 Lac, making some assump- 
tions and restrictions not supported by our study. In the 
present paper, we discuss our conclusions w hich differ from 
those of IDziembowski fc Pamvatnvkhl (|2008l ) . 



2 OBSERVATIONS AND DATA REDUCTION 

The data originate from a spectroscopic multisite campaign 
dedicated to o ur studied star . In a ddition, we added the 
903 spectra of iMathias et al.l (| 19941 ) , which allowed us to 
increeise the time span and, thus, to achieve a better fre- 
quency accuracy. In total 1820 observations were gathered 
using 9 different small- to medium-sized telescopes spread 
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Figure 2. The Sim Une profiles of 1 night (JD 2452898) taken 
with the McDonald telescope. The arrow indicates the fiux scale 
of one spectrum. 



over the northern continents. Table [T] summarizes the log- 
book of our spectroscopic data. The resolution (A/ A A) of 
the instruments ranged from 30 000 to 80 000 and the aver- 
age S/N ratio near 4500 A between 180 and 380. 

All data were subjected to the normal reduction pro- 
cess, which consists of de-biasing, background subtraction, 
flat-fielding and wavelength calibration . Finally, the helio- 
centric corrections were computed, and all spectra were nor- 
malized to the continuum by fitting a cubic spline function. 

For our study, we considered the Sim triplet around 
4567 A because its characteristics simplify the modelling 
of the line-profile variations, which we use for mode iden- 
tification. Indeed, these silicon lines are sufficiently strong 
without being much affected by blending. Moreover, they 
are dominated by temperature broadening, such that the 
intrinsic profile can simply be modelled with a gaussian. Fi- 
nall y, they are not too sensitive to temperature v ariations 
(see iDe Ridder et al.|[200^ : lAerts fc De Catll2003l ). so that 
neglecting them remains justified. 

The similar average radial velocity of about —15 kms^^ 
for 1990-199 2, 2003-2004 and earlier measurements (e.g. 
IWilsonI Il953t ) definitely excludes the possibility of 12 Lac 
being a spectroscopic binary. Besides the pulsation effect, 
the different centroids are due to the different zero points of 
the different telescopes. Before extracting the pulsation in- 
formation from the whole dataset, it is necessary to correct 
for this effect. The spectra were shifted in such a way that 
the constant of a least-squares sine fit using the first three 
dominant modes is put to the same value for each observa- 
tory. 

In Fig.[T]we superimpose the average profiles of the Si ill 
A4553 A line computed for each observatory separately, af- 
ter correction for the different zero points. We note that the 
centroid of the lines are in good agreement but not exactly 
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Figure 3. The densest part of the time series of the Sim radial velocities ((f^)) of 12 Lac derived from each spectrum taken during the 
dedicated multisite campaign (7 Aug. 2003 - 17 Nov. 2003). 



Table 1. Log of our spectroscopic multisite campaign. The Julian Dates are given in days, AT denotes the time-span expressed in days, 
A'^ is the number of spectra and S/N denotes the average signal-to- noise ratio for each observatory measured at the continuum between 
4500 and 4551 A. 



Observatory 


Long. 


Lat. 


Telescope 


Julian Date 


Data amount and quality 


Observer(s) 


(Name of the instrument; resolution) 








Begin 


End 


AT 


N 


S/N 




2003-2004 








-2450000 










Apache Point Observatory 


-105°49' 


-|-32°46' 


3.5-m 


2927 


2928 


2 


93 


300 


JK 


(ARC; 31500) 




















Bohyunsan Astronomical Observatory 


-l-128°58' 


-l-36°09' 


1.9-m 


3109 


3109 


1 


27 


180 


DM 


(BOES; 30000) 




















Dominion Astrophysical Observatory 


-123°25' 


-|-48°31' 


1.2-m 


2858 


2859 


2 


66 


250 


SY 


(45000) 








2909 


2911 


3 


42 






McDonald Observatory 


-104°01' 


-1-30° 40' 


2.1-m 


2893 


2901 


9 


210 


380 


GH 


(Sandiford; 60000) 






2.7-m 


3195 


3198 


4 


31 


350 


GH 


Nordic Optical Telescope 


-17° 53' 


+28°45' 


2.6-m 


2947 


2951 


5 


32 


370 


KU, JT, II 


(SOFIN; 80000) 








2957 


2960 


4 


47 






Okayama Astrophysical Observatory 


-fl33°35' 


-1-34° 34' 


1.88-m 


3201 


3202 


2 


7 


180 


SM 


(HIDES; 685000) 








3205 


3205 


1 


6 














3208 


3210 


3 


17 






Thiiringer Landessternwarte Tautenburg 


-|-11°42' 


-1-50° 58' 


2-m 


2863 


2869 


7 


103 


240 


HL 


(67000) 








2886 


2889 


4 


100 














2895 


2900 


6 


105 














2914 


2918 


5 


31 






1990-1992 








-2440000 










Observatoire de Haute-Provence 


+5°42' 


-l-43°55' 


1.52-m 


8135 


8135 


1 


70 


220 


PM 


(AURELIE; 42000) 








8224 


8224 


1 


52 














8227 


8227 


1 


35 














8845 


8858 


14 


582 














8884 


8889 


6 


164 






Total 












82 


1820 







at the same position. Larger differences can be seen for the 
width and the depth of the Unes. Such deviant average pro- 
files are due to the limited time spread of the observations 
together with the multiperiodic character of the pulsations. 
In the case of a simil ar spectrosc o pic m ultisite campaign 
devoted to v Eridani, lAerts et al.l l|2004l ) observed similar 
differences, which are not caused by instrumental and/or 
signal-to-noise ratio effects. 

Fig. [2] displays the large and complex radial-velocity 
variations of 12 Lac. The plotted spectra were taken with 



the McDonald telescope during only one night. In Fig.O we 
show the radial velocities for the densest part of the multi- 
site campaign (7 Aug. 2003 - 17 Nov. 2003). A clear beating 
pattern is seen. 

3 FREQUENCY ANALYSIS 

We performed a frequency analysis on the first three ve- 
locity moments {v^), (w^) and (li'^) (see lAerts et al.lll992l . 
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Table 2. Frequencies and radial velocity amplitudes of the 
first moment of the Si III A4553 A line together with their S /N 
ratio (we refer to the text for e xplanation). Error estimates 
llMontgomerv fc 0'Donoghuelll999i ) for the independent frequen- 
cies range from ±0.000002 d"^ for /i to ±0.00002 d^^ for /g. The 
error on the amplitude is 0.01 kms~^. 



ID 



Frequency 







[d-1 ] 


[mHz] 


[kms ^ ] 




/l 




5.178964 


59.941713 


14.50 


99.6 


h 




5.334224 


61.738704 


7.70 


52.2 


h 




5.066316 


58.637917 


6.26 


42.7 


k 




5.490133 


63.543206 


2.61 


17.5 


fa 




0.342841 


3.968067 


1.34 


7.1 


h 




4.241787 


49.09475 


0.91 


6.7 


h 




5.218075 


60.394387 


0.84 


5.8 


fr 




6.702318 


77.573125 


0.62 


4.4 


fs 




7.407162 


85.731042 


0.67 


4.9 


h 




5.84511 


67.65184 


0.79 


5.2 


2/1 




10.35814 


119.88590 


0.63 


4.5 


/2 4 


/3 


10.40056 


120.37693 


0.72 


3.8 


/i4 


/2 


10.51319 


121.68044 


2.59 


19.7 


2/i 


+ /3 


15.42400 


178.51860 


0.72 


6.6 


/i4 


/2+/3 


15.57950 


180.31838 


0.71 


6.4 



for a definition of the moments of a line profile) of 
the Si m A4553 A line by means of the program Pe- 
riod 04 (iLenz fc Bregej |2005|). For some 13 Cephei stars 
fseelTelting et aLlllQQTt ISchriivers et al.ll2004l : iBriquet et al.1 
I2OO5I ) a two-dimensional frequency analysis on the spec- 
tral lines led to additional frequencies compared to a one- 
dimensional frequency search in the moments. For 12 Lac 



ID analysis on (w"*) and {v^) did not reveal additional inde- 
pendent periodicities compared to those present in the first 
moment. We also do not find additional frequencies in the 
2D frequency analysis. The 2D analysis is more sensitive to 
the detection of high-degree modes than a ID analysis so 
here we have already an indication that the frequencies are 
low-degree modes. 

In what follows we only describe our frequency analysis 
on (w^), which is the radial velocity from which the fitted 
averaged radial velocity is subtracted. 

The frequencies were determined with the standard 
method of prewhitening the data. In addition we made use 
of a procedure available in Period04 which allows to improve 
the detected frequency from non-linear least-squares fitting 
with the maxima in the Fourier transform as starting val- 
ues. At each step of prewhitening we subtracted a theoretical 
multi-sine fit with the amplitudes, phases and also optimized 
frequencies that yielded the smallest residual variance. All 
peaks exceeding an amplitude signal -to-noise r a tio ab ove 4 
in the Fourier periodogram (,Brcacr et al.l Il993l . Il999l ) were 
retained. The noise level was calculated from the average 
amplitude, computed from the residuals, in a 2d~^ interval 
centered on the frequency of interest. 

The detected frequencies are listed in Table [2] and the 
Fourier periodograms are shown in Fig.|4] Our frequency 
analysis reaffirmes the presence of the already well-known 
5 main frequencies (|Mathias et al.|[l994l ) together with the 
new indepe ndent signals a nd co mbination frequencies dis- 
covered by [Handler et al.l (|2006l ) . We note that one addi- 



Amplitudc in {v^) S/N [ 



tional independent frequency [fp — 5.30912 d~^) was found 
by the latter authors and not by us. An important result 
is that we clearly recover the low-frequency signal possibly 
originating from a g-mode. Such a SPB-like oscillation does 
not seem to be uncommon among /3 Cephei stars. For in - 
stance, it was also observed i n 19 Mon (iBalona et al.ll2002l) 
in V Eridani (|Handler et al.l |2004| ). in S Ceti (|Aerts et all 
200^ and in V1449 Aql (Briquet et al., in preparation). We 
also point out that such hybrid SPB-/3 Cep pulsators are the- 
oretically predicted for ho tter stars than prev iously thought 
jMiglio et al.ll2007l versus |Pamvatnvkh|[l999l ) . 

The equivalent width (EW) of the silicon line of 12 
Lac clearly varies with the two dominant modes /i and /2. 
The sum frequency /i + /2 is even present. Fig. [5] represents 
phase diagrams for both frequencies. Such a strong EW vari- 
ability of 9% is generally not observed in pulsating B-type 
stars. Indeed, typical relative EW p eak-to-peak amplitudes 
are of the order of a few percent (|De Ridder et all I2OO2I : 
[d7 Cat fc Aertsll20()^ ). However, for several (3 Cephei stars, 
one clear sinusoidal variation is observed in their EW, with 
the same frequency as in the radial velocity. The known 
cases are targets pulsating multiperiodically with a high- 
ampl itude domin ant radial mode, s uch as & Ceti (|Aerts et al.l 
Il992h . V Eridani (|Aerts et al.ll2004h and V1449 Aql (Briquet 
et al., in preparation). In the case of 12 Lac, we note that the 
EW signal is multiperiodic with a dominant (1, l)-mode, fol- 
lowed by a radial mode of lower amplitude (see next section) . 
For the four strongest modes, the frequency values derived 
from our spectroscopic measurements are ex actly the same, 
within the errors, as the photometric ones bv [Handler et all 
(|2006l ). For the modes with radial velocity amplitudes typ- 
ically lower than 1 kms~^, there are differences between 
values derived from both kinds of data, which might have 
some infiuence on the mode identification outcome. Because 
the number of observations as well as the time span of the 
datasets in the literature are much larger than the dedi- 
cated spectroscopic ones, the literature frequency values are 
expected to be more accurate and hence we adopt them. 



4 MODE IDENTIFICATION 
4.1 The methods used 

We used two different and independent me thods to identify 
the m odes, namely the moment method (|Briquet fc Aertj 
|2003| ) a nd of the Fourier parameter fit method (FPF 
method. [Zima|[20061 . 

The basic principles of both methods are the follow- 
ing. With the moment method the wavenumbers (l,m) and 
the other continuous velocity parameters are determined in 
such a way that the theoretically computed first three mo- 
ment variations of a line profile best fit the observed ones. 
In the FPF method, for every detected pulsation frequency, 
the zero point, amplitude, and phase are computed for ev- 
ery wavelength bin across the profile by a multi-periodic 
least-square fit. Subsequently the zero-point profile, the am- 
plitude and the phase values across the profile are fitted with 
theoretical values derived from synthetic line profiles. The 
first technique has the advantage to be less computation- 
ally demanding, allowing to test a huge grid for [t, m) and 
for the other parameters. Moreover, the modes are deter- 
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Figure 4. Amplitude spectra of 12 Lac computed from tlic first moment of the Si III A4553 A lino. The uppermost panel shows the 
spectral window of the data. All subsequent panels show the periodograms in different stages of prewhitening. The significance limit in 
red is calculated according to 4 times the noise over a 2d~^-range around the significant frequency. 



mined simultaneously by fitting a multiperiodic signal tak- 
ing into account the coupling terms appearing in the second 
and third moments. For the second technique, the fitting is 
carried out by applying genetic optimization routines in a 
large parameter space and mono-mode fits are used in order 
to speed up the computations. Once the best solutions are 
constrained, a multi-mode fit is performed. The strength of 
the FPF method is its ability to estimate the significance 
of the derived mode parameters by means of a te st (all 
y^-va lues listed in the text are reduced values, see IZimal 
l2008l ). 

The FPF method is not optimally suited to study high- 
amplitude modes, but as described further, the applica- 
tion of the F PF method f or 12 Lac pro ved to be suc cesful. 
As shown bv IZima et al] (|2004l ') and by IZimal (|2006l '). both 
mode identification techniques constrain the value of the az- 
imuthal number m better than the degree £ for low vsini. 
Since the photometric data pr ovided us with unam biguous 
1- values for the 5 main modes (|Handler et al.|[2006l ). a valid 
strategy is to make use of spectroscopy to identify the val- 
ues of m, adopting the 1-values obtained from photometry. 



Such a methodology a lready proved to be successful for two 
other f3 Cephei stars (|Briauet et al.|[2005l : iMazumdar et al.l 
I2OO6I ). Mode identification results have been obtained with 
the software package FAMIA^ (,Zima 2008l ). 

In what follows we present in detail our mode iden- 
tification by the FPF method. The moment method gave 
compatible results, but did not give additional constraints 
on the {£, m) than the ones we present here. Both techniques 
found the same m-values for the three main modes, giving 
us much confidence in our outcome. For the fourth mode, 
the moment method derived the sign of m in agreement with 
the FPF method, which could further constrain its value. 

4.2 Derivation of the wavenumbers 

For the FPF method, a sophisticated theoretical formalism 
for the modelling of the displacement field is used and in- 

^ FAMIAS has been developed in the framework of the FP6 Euro- 
pean Coordination Action HELAS - |http : //www.helas-eu. org/ [ 
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Figure 5. Phase diagrams of tlie equivalent widtli of the Sim 
A4553 A line, for /i (top panel), and for /2 (bottom panel) after 
prewhitening with /i . 



eludes, for instance, first-order effects of the Coriolis force 
as well as temperature variations of the stellar atmosphere. 
We refer to lZima (,2006, . ,200^ for details on the computation 
of the synthetic line-profile variations. The fixed parameters 
are the mass and the radius. When taking into account tem- 
perature effects, the effective temperature Toff and the loga- 
rithm of the gravity log g have also to be fixed. For 12 Lac , 
we adopted t he values M = 13.5 ([Handler et al.ll2006t ). 
R = 6.9 J^a dHaiidler et al.ll2006l ). T^a = 24500 ± lOOOK 
(|Morel et al.ll2006h . log o = 3.65 ± 0.15 (|Morel et al.ll2006h . 
We note that the mode identification results are robust when 
using different values for these parameters, within the errors. 

In a first step we identified the 3 main modes (/i,3,4) by 
fitting mono-mode Fourier parameters, which are computed 
from synthetic line profiles evenly sampled over one pulsa- 
tion cycle (for more explanations see[Zima 20061). Moreover, 
we neglected the equivalent width variations of the spec- 
tral line due to temperature variations. In order to restrict 
the parameter space, we estimate roughly the values of the 
projected rotational velocity vsmi, the width of the intrin- 
sic gaussian profile cr, and the equivalent width W from a 
least-squares fit of a rotationally broadened synthetic pro- 
file to the zero-point profile. This gives only a crude estimate 
of these three parameters, because the pulsational broaden- 
ing is ignored in this fit, but it serves the purpose of defin- 




-120-100 -80 -60 -40 -20 20 
Velocity [km s~^] 

Figure 6. Observed zero-point profile (points) and best fit when 
taking into account the four main modes (solid line) by means of 
a synthetic rotationally broadened profile (iisini = 36.7kms~^, 
W=15.3kms-\ a = 22.1 kms'^). 



ing our appropriate grid. Moreover v sin i will be overesti- 
mated. When taken into account the pulsational broaden- 
ing caused by the four main modes we acquire a vsmi of 
36.7 km s^^ (see Fig.|6]). A genetic optimization with the fol- 
lowing free parameters was afterwards adopted: I £ [0, 4] 
with a step of 1, m G \—l,t\ with a step of 1, the sur- 
face velocity amplitude a £ [10, 110] kms"'^ with a step of 
1 kms"'^, the stellar inclination angle i £ [5, 85] ° with a step 
of 5 °, V sin i £ [20, 48] kms~^ with a step of 0.5 kms~^ and 
<J G [10, 30] kms~^ with a step of 0.1 kms~^. We fixed /2 as 
a radial mode and only fitted the continuous parameters for 
that mode. 

In Table [S] the 5 best solutions from a fit to the 
Fourier parameters are listed for each frequency. The 95%- 
confidence limit of is 1.15. The most probable solutions 
are (/i,mi) = (1,1), {h,m3) = (1,0) and (l4,,m4,) = (2,1) 
(where a positive m- value denotes a prograde mode). We 
point out that the agreement between theory and observa- 
tions is not perfect. More precisely, the observed asymmetry 
with respect to the centroid of the amplitudes across the line 
of the dominant modes is not completely reproduced by our 
model. Such an asymmetry can be caused by EW varia- 
tions of th e intrinsic line profile due to local temperature 
variations (|Schriivers fc Teltind [l999l ) . We considered EW 
variations in our fitting procedure but we did not obtain an 
improved solution compared to the case of a constant EW. 

In a second step we attempted to improve the solutions 
by simultaneously fitting the four main modes and using all 
the original data points. This however was computationally 
not possible. To restrict the computation time we took a sub- 
set of our complete time series. Our subset resulted in 388 
spectra with a time span of 4 months. We also fixed the wave 
numbers according to our mono-mode fit described above. 
We succeeded in improving the surface velocity amplitudes 
together with the inclination, usini and the width of the 
intrinsic profile. The values are lower than in the mono- 
mode fit and we are able to reproduce part of the asymmetry 
in the amplitude across the line as explained above (see Ta- 
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Table 3. Mode parameters derived from the FPF mono-mode method. For each pulsation mode, the five best solutions are shown together 
with the second best {£,m) identifications. For /s we list the different solutions for m. as is the surface velocity amplitude {kms~^); i is 
the stellar inclination angle in degrees; v sin i is the projected rotational velocity, a is the width of the intrinsic profile, both expressed in 
kms~^. Due to spherical symmetry we have not indicated the inclination angle for radial modes. The 95% significance limit is at = 1-15. 
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Table 4. Mode parameters derived from the FPF method 
through a simultaneous fit of the four main modes. We used the 
same symbol conventions as in Table|3] The 95% significance limit 
is at = 1-31. 
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ble|4]). Fig. [7] illustrates the best models for the amplitude 
and phase across the line profile for the four main modes. 

Finally, we tried to identify the other low-amplitude 
modes with the same methodology and by imposing the 
wavenumbers {£, m) deduced for the 4 main modes. Unfor- 
tunately, we could not obtain additional conclusions. Our 
failure is not very surprising since these modes have am- 
plitudes lower or of the same order as the harmonics and 
combination frequencies also present in the signal while non- 
linear terms are not taken into account in our line-profile 
modelling. Fig.[S]depicts the amplitude and phase across the 
line for these low amplitude modes. Despite the fact that the 
phase across the profile is well determined for these modes 
and seems to point to axisymmetric or prograde modes, the 
values between the different options were not discrimi- 
native so we do not use that information in the following. 
In Table [SI we summarize the fi nal outcome of both the 
photometric (|Handler et al.1 [200^ and spectroscopic mode 
identifications. 



Table 5. Final results for the mode identifications of 12 Lac 
from our spectroscopic analysis t ogether with the re sults from 
the photometric amplitude ratios l lHandler et al.ll2006ft . 
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1,2 






5.30912 
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4.3 Derivation of the surface equatorial rotational 
velocity and the inclination 

From mode identification we can also estimate the inclina- 
tion, V sin i and the equatorial rotational velocity Wcq. By us- 
ing Bs, a weight we can construct histograms for i, vsvai 
and «cq. In Fig.|9]we display these histograms. They are com- 
puted by considering all the solutions of the FPF method 
with the correct wavenumbers of the four dominant modes 
through a simultaneous fit, and by giving each parameter 
(ifc, vsmik and Heq.fe) its appropriate weight Wk = Xo/xl, 
where Xo is the x^-value for the best solution. By calculating 
a weighted mean and standard deviation of the data, we get 
i = 48±2°, vsini = 36±2kms"^ and Ueq = 49±3kms~^ 
This leads to a surface rotational frequency between 0.11 and 
0.12 d^^ for the appropriate stellar models in Table[6l which 
will be discussed further on. In the same way we constructed 
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Figure 7. Observed Fourier parameters (black points with error bars) across the line profile together with the best fit for the four 
identified modes (full line). The amplitudes are expressed in units of continuum and the phases in 2-k radians. 
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Figure 8. Observed Fourier parameters (black points witli error bars) across the line for the detected pulsation frequencies (fgjy-g). 
The amplitudes are expressed in units of continuum and the phase in 2it radians. 



these histograms with the solutions from the multi-mode fit 
of the moment method. However, the moment method did 
not have the ability to limit the range for i and, by implica- 
tion, leads to a less trustworthy estimate of Vcq. The reason 
for this is probably that the moment method is based on 
integrated quantities over the line profile, while in the FPF 
method we use all the information across the line profile 
as a whole, which is more sensitive to the position of the 
nodal lines across the stellar surface and thus to the stellar 
inclination. 



vective overshooting parameter. Such asteroseismic con- 
straints have recently been deri ved for a fe w /3 Cephei 
stars : 16 Lac (iThoul et al.1 12003| ) , V836 Cen (lAerts et all 
20031: iDupret et all 120041') . v Eri (IPamvatnvkh et al. |2004| 



Ausseloos et al.l 



gooj r [Pziembowski fc Pamvatnvkhl 



5 Ceti ( Aerts et aLlbOOel). g CMa (iMazumdar et a" 



e Oph (IBriauet et al.ll2007^. Finally, I AusseloosI i2Qok) and. 



200S), 



20061 ). 



more recently, iDziembowski fc Pamvatnvkh "(|2008l ) already 
made a comparison with stellar models for 12 Lac. 



5 MODELLING 

In this section, we check if state-of-the-art stellar mod- 
els with standard physics can account for the observed 
frequency spectrum together with the derived wavenum- 
bers {£, m) for our studied star. By doing so, we can 
constrain 12 Lac's model parameters, which are the 
mass, the central hydrogen abundance and the core con- 



5.1 Numerical tools and physical inputs 

In our work, we used the following numerical tools and 
physical inputs. The stellar models for non-rotating stars 
were computed with the evol utionary code CLES (Code 
Liegeois d'Evolution Stellaire, IScuflaire etal] l2008bh. We 
used the QPAL2001 equation o f state ( Rogers fc Navfonovl 
l2002l : ICaughlan fc Fowldll988t) . with nuclear reaction rates 
from Formicola et al.l (|2007|) for the "N(p,7)^'^0 cross- 
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Figure 9. Histograms for the inclination, projected rotation ve- 
locity and equatorial rotational velocity of the star derived from 
the FPF method. 



section. Convective transport is treated by usinR the clas- 
sical Mixing Length Theory of convection (|Bohm-Vitensel 
1 19581 ). Fo r the chemical compos ition, we used the solar mix- 
ture fromlAsplund et all lj2005l l . We used OP opacity tables 
ISeatonlbOOsF computed for this latter mixture. These ta- 
bles are completed at logT < 4 .1 with the low temperature 
tables of iFerguson et all (|2005h . 

Our choice for the metal mixture and t he opacities 
is just ified by the following arguments. First, I Morel et al.l 
l|2006l ) found that the abundances of He, C, N, O, Mg, Al, Si, 
S and Fe are, within the errors, in good agreement with the 
solar values derived from ti me-dependent , three-d imensional 
hydrodynamical models of lAsplund et all (|2005l . and refer- 
ences therein) . Secondly, lAusselood (|2005l ) performed a pre- 
liminary seismic study of 12 Lac and was confronted with the 
problem that, when using OPAL opacities, not all observed 
pulsation modes can be excited for realistic values of the 
metallicity and core overshooting p arameter. This excita - 
tion problem is en countered with t helCrevesse et al.l (|l993l ) 
but also with the lAsplund et al] (|2005h mixture, even for 
large values of Z (typically larger than 0.016). If we con- 
sider values of Z smaller than 0.016, which are more rep- 
resentative of t he star, the problem is obviously increased. 
As discussed in lBriguet et al.l (|2007) for 9 Oph, modes that 
are stable with OPAL o pacities ( Rogers fc Iglesia^ 1 19921 : 
[Rogers fc Navfonovl [20o3 ) can be predicted to be unstable 



with OP opacities. The reason is that in a typical /3 Cephei 
model the OP opacity is 25% larger than the one for OPAL 
in the region where the driving of pulsation modes occurs. 
We consequently wish to check if the use of OP opacity ta- 
bles improves t he situation in the ca se of the modelling of 12 
Lac. We refer to lMigho et"al] (|2007l ') for a detailed discussion 
on the implication of the adopted opacity tables and metal 
mixtures on the excitation of pulsation modes in the case of 
B-type stars. 

For each stellar model, we calculated the theoretical fre- 
quency spectrum of low-order p- and gr-modes with a degree 
of the oscillation up to ^ = 4 usin g a standard adiabati c 
code for non-rotating stellar models jScuflaire et al.ir2008al ). 
which is much faster than a non-adiabatic code but leads 
to the same theoretical pulsation frequencies within the 
adopted precision of the fit, which was 10"^ d~^ Once the 
models fitting the observed modes are selected, we checked 
the excitation of the pulsation mo des with the l inear non- 
adiabatic code MAD developed bv iDupretl (|20o3). 

5.2 Seismic analysis 

Besides the pulsation characteristics summarized in Ta- 
bleO we also use Tea, logg and the metallicity Z 
as additional observational constraints. Geneva photom- 
etry of 12 Lac provides Teg = 23500±700 K and 
log .9 = 3.4±0.4 jHandler et al.ll2006l ) . From lUE spectra, 
iNiemczura fc Daszviiska-Daszkiewicj (|2005|) derived T^s = 
23600±1100 K and logg = 3.65 ± 0.15. iMorel et al.1 (120061 ) 
obtained T^s = 24 500±1000 K and logg = 3.65±0.1 from 
optical spectra. Us i ng an o ther method on the same spec- 
tra as iMorel et al] (|2006h , iLefeveij (|2008|) deduced T eff = 
24000±1000 K and log g = 3.60±0. 1. [Morel et al.1 l|2006l ) also 
determined a metallicity Z = 0.0089±0.0018. 

Stellar models are parametrized by the initial hydrogen 
abundance X, the core convective overshooting parameter 
Oov, the metallicity Z and the mass M. We adopted X = 
0.72 but using a different value for X in a reasonable range 
does not change the conclusions. 

We considered two values for Z, i.e. 0.010 and 0.015. 
The first chosen value corresponds to a metallicity typical 
of B dwarfs in the solar neigh bourhood, also fou nd for 12 Lac 
within the errors (see above, iMorel et al.ll200^ ) . The second 
choice corresponds to the solar composition and the reason 
for also adopting this value is the following. The star being a 
relatively slow rotator and its surface convective zone being 
very thin, diffusion mechanisms can alter the photospheric 
abundances, leading to Z smaller at the surface compared 
to the interior of the star. However, since these effects are 
limited to the very surperficial layers, including a diffusion 
treatment or not does not affect the derived stellar parame- 
ters. It has been carefu l ly che cked in the case of 6 Oph. We 
refer to iBriguet et al] (|2007l ) for a deeper discussion. For 
simplicity, we computed non-rotating stellar models with- 
out taking into account diffusion, but we also considered a 
Z higher than what is derived from spectra. 

Several values for Oov in the range [0.0-0.5] were tested. 
The reason is that recent asteroseismic modelling of several 
P Cephei stars (listed above) showed the occurrence of core 
overshooting with Oov in the range [0.05-0.25] for all studied 
stars except 9 Oph for which an even higher value of ^ 0.45 
was determined. Moreover, other works corroborate the ne- 
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cessity to include core over shooting in mod elling of massive 
B-type stars. For instance, iDeuprej (|2000l ) derived a value 
of about 0.45 by means of 2D hydrodynamic simulations of 
zero-age main-sequence convective cores. We also mention 
that, wit h his s t udy o f 13 detached double-lined eclipsing 
binaries, IClared |20o3) found that models with core over- 
shooting of aov ~ 0.2 are needed to analyse ~10 Mq stars. 

The observed radial mode with /2 = 5.334224 d~^ and 
the i = 1 zonal mode with /s = 5.066316 d~^ were selected 
as the first two frequencies to be confronted with the theo- 
retically predicted frequencies. This was an obvious choice 
since these are zonal modes (m = 0) and thus we do not need 
any information about the equatorial frequency of rotation 
to fit them with the model frequencies. We considered two 
different scenarios: /2 being either the radial fundamental 
mode or the first overtone. 

The theoretical frequency spectrum with /2 being the 
radial fundamental mode is globally compatible with the ob- 
served frequencies of 12 Lac. However, all such models fail 
to reproduce the £ = 1 mode with /y = 6.702318 d^^. If /2 
corresponds to a radial fundamental mode, /r can only be 
an ^ = 2 mode according to theoretical models, while the 
observed phot ometric amplitude r atios for this mode point 
out an £ = 1 (|Handler et al.ll2006h . Unfortunately, we can- 
not confirm this strong constraint since our spectroscopic 
mode identification was not conclusive for this mode. How- 
ever, stellar models with /2 corresponding to the radial first 
overtone can account for a, £ = 1 mode with /y. We conse- 
quently favour a radial order n = 2 for /2. An additional 
strong argument to favour /2 as the radial first overtone 
is the following. The Ledoux rotational splitting is defined 
by /Ledoux = rn0nif rot, whcrc Pni is a structure constant 
depending on the stellar model, and originates from the fre- 
quency splitting defined as 

U = fo + mP^iLot =yo + m r K{r)^%, (1) 

Jo ^'^ 

where /m is the cyclic frequency of a mode of azimuthal- 
order m, 0,{r) is the angular velocity, and K(r) is the ro- 
tation al kernel (for more information, see IScuflaire et al.l 
l2008al ). With /2 as fundamental, the difference between the 
frequency /s = 4.241787 d^^ and the nearest quadrupole 
mode (at ~3.8d~^) is too large to be explained by the 
Ledoux rotational splitting and thus /s cannot be fitted. 
With /2 as first overtone, /s is easily identified as the g2 
mode with {£,m) = (2, 1) as shown below (see also Fig. lllf) . 
With /2 being the first overtone, /s is identified as pi. 

The couple {X, Z) being fixed, the fitting of two fre- 
quencies suffices to derive the stellar age and stellar mass 
for adopted Qqv. The positions in the log(TefF) — logg dia- 
gram of the models matching /2 and /a for different values 
of Oov are shown in Fig llOl The corresponding parameters 
are listed in Table[6]for Z = 0.015. The mass decreases when 
the core overshooting parameter increases. We also note the 
perfect agreement between the \ogg of the models and the 
value deduced from spectroscopic observations. Moreover, 
the best agreement in Tea corresponds to models with Qov 
between 0.0 and 0.4 and a mass between 10.0 and 14.4 Mq. 
Note that a decrease in Z implies an increase in mass and, 
considering a Z of 0.010 instead of 0.015 increases the mass 
by 0.5 Mq at most. 



Table 6. Physical parameters of the models that match /2 (being 
the first overtone) and /a, with X = 0.72 and Z = 0.015. Xc is 
the central hydrogen abundance. The age t is expressed in million 
years. 
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0.1 


13.1 


24500 


3.68 
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Figure 10. The error box represents the position of 12 Lac in 
the log(T(,fl ) — logg diagram. The positions of the models which 
fit exactly /2 (being the first overtone) and /s are also shown for 
different «ov values. 



5.3 Radial order identifications 

The frequency /i = 5.178964 d~^ can only belong to the pi 
triplet. The frequency splitting /i — fs is 0.1126 d^^. Ac- 
cording to the models, = 0.5 for this mode. This yields 
an averaged frequency of rotation /rot = (/i — f3)/m/3ni- 
This averaged rotation frequency leads to a rotational ve- 
locity twice as high as the surface value estimated from our 
spectroscopic observations (with corresponding velocities of 
lOOkms"^ versus 49kms~^). We thus confirm that 12 Lac 
must rotate more rapidly in it s inner parts than at the sur- 
face, a s already emphasized bv lDziembowski fc Pamvatnvkhl 
(|2008l ). 

As can be seen in Fig. 1111 the {£ = 1,P3) mode can 
account for /y = 6.702318 d^^. The match is not perfect for 
all models in Table|6] This is probably due to the fact that 
there can occur small shifts due to rotation. On the other 
hand a change in metallicity will also cause a small shift 
of the theoretical frequencies (see Sect. 15. 5| ). The frequency 
fs = 7.407162 d-^ and = 5.84511 d"^ can only be the 
{£ = 2,p2) mode and the {£ = 1,P2) mode, respectively. 
Indeed, no theoretical frequency of an ^ = 1 zonal mode 
and an £ = 2 zonal mode is in the vicinity of /g and /g, 
respectively. We also note that an / = radial-mode solution 
for fs is excluded (see Table[5l). 

The frequency fp = 5.30912 d~^ found in photometry 
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Table 7. The I, n and m identifications according to our mode 
identification and modelling with the radial mode as first over- 
tone. We use the usual convention that negative n- values denote 
g-modes (see e.g. lUnno et al.l ll989l). 



Table 8. Observed photometric amplitude ratios taken from 
[Handler et"al] | |2006D . We also list our best fitting theoretical am- 
plitude ratios. They originate from the model with Qov = 0.4 
given in Table|6] 



ID 


Frequency 
[d-i ] 


I 


n 


m 


/i 


5.178964 


1 


1 


1 


h 


5.334224 





2 





h 


5.066316 


1 


1 





U 


5.490133 


2 





1 


h 


4.241787 


2 


-2 


1 


h 


5.218075 


4 


-2 or -1 




fr 


6.702318 


1 


3 




fa 


7.407162 


2 


2 




h 


5.84511 


1 


2 






5.30912 


2 









ID 


Frequency 


v/u 




V/u 






[d-i ] 


obs. 


theory 


obs. 


theory 


/2 


5.334224 


0.529(7) 


0.64 


0.456(7) 


0.56 


/3 


5.066316 


0.716(8) 


0.73 


0.686(8) 


0.67 


/4 


5.490133 


0.82(1) 


0.83 


0.78(1) 


0.78 



and the frequency fi — 5.490133 were both identified as 
quadrupole modes and can be accounted for by the {I — 2, /) 
mode. Indeed, the theoretical Ledoux frequency splitting as- 
suming the averaged internal rotation frequency value /rot 
derived from /i — /a ranges in [0.186-0.190] d~^ for differ- 
ent overshooting values, which corresponds to the difference 
between fp and /i. The frequency fp thus corresponds to 
{£, m) = (2, 0) and f^ to (I, m) = (2, 1) in accordance with 
the m— identification for fi. Note that the theoretical fre- 
quency value for the [l — 2, /) mode is in best agreement 
for Qov = 0.0 and is 5.287 d~^ 

/s — 4.241787 d~^ being a prograde £ — 2 mode, is 
uniquely identified as the {£ = 2,g2) mode. The corre- 
sponding theoretical zonal mode has a frequency between 
4.048 and 4.070 d~^ with a Ledoux splitting for that mode 
ranging between 0.183 and 0.185 d~^ deduced from /rot and 
Pni = 0.80. In addition, we thus conclude that m = 1 for 
this mode. 

/e — 5.218075 d"'^ may be either identified as {£ = 
4, g2) or {£ — 4, gi) (see also Fig. llip . The models give 
Pni = 0.97 for this mode. In the first case, the zonal fre- 
quency ranges between 4.38 and 4.55 d~^ with a Ledoux 
splitting of 0.21-0.22 d"^ deduced from frot- This corre- 
sponds to an m-value of 3 or 4. In the second case, the zonal 
frequency ranges between 5.57 and 5.70 d~^ with a Ledoux 
splitting of 0.20-0.21 d^^. This corresponds to an m-value 
of -2. 

Finally, nothing can be concluded for fg = 0.342841 d^^ 
because of the dense theoretical frequency spectrum in the 
low-frequency range and because of the lack of observed con- 
straints on the wavenumbers of this mode. We additionally 
note that we cannot definitely exclude this frequency to be 
linked with the frequency of rotation. For instance, /g/3 
is compatible with our obs erved surface equatorial rot ation 
frequency. As explained in iBriguet et al. 1 (|200ll . l2004l ). one 
way to discriminate between the pulsation and rotation in- 
terpretation is to compare the variability of the considered 
frequency in lines of different chemical elements. Unfortu- 
nately, we could not achieve such a test because lines other 
than Si lines are at our disposal only for two nights of ob- 
servation. This time span is not long enough to recover such 
a low frequency value. 



5.4 Comparison with photometry 

We used our best models listed in Table[5] to compare 
our results w i th th e photometric mode identification of 
[Handler et al.l (|2006t ). We computed theoretical photomet- 
ric amplitude ratios for our three main identi fied pulsation 
modes (/2, /a, /4; /i belongs to /a) following iDupret et al.l 
(|2003l ). We computed the amplitude ratios for Stromgren it 
and V, and for Johnson V. For this computation we must de- 
termine the nonadiabatic pulsation mode parameters fr and 
fg, where /t corresponds to the local effective temperature 
variation and fg corresponds to the local effective gravity 
variation, both at the level of the photos phere. We used the 
non-adiabatic code MAD (|Dupretll200ll ) to compute these 
fr.g parameters for our best seismic models of 12 Lac. 

We list the observed and theoretical photometric am- 
plitudes for the three main m odes in Table[8l We refer to 
Fig. 4 of [Handler et al.l (|2006l ') for the results of the photo- 
metric mode identification. For all the models, we observed 
that the amplitude ratios decrease with decreasing mass. For 
/2 we find amplitude ratios for v/u ranging from 0.68 until 

0. 64, and for V/u ranging from 0.51 u ntil 0.56. Thus , we en - 
counter exactly the same problem as [Handler et al] (|2006h . 

1. e., we could not find an perfect agreement for the observed 
ratios of /2. For /a we find amplitude ratios for v/u ranging 
from 0.84 until 0.73, and for V/u ranging from 0.74 until 
0.67. Our models give amplitude ratios for /4 for v/u rang- 
ing from 0.90 until 0.83, and for V/u ranging from 0.84 until 
0.78. Our conclusion is that we can reproduce the observed 
photometric amplitude ratios for the two main modes (/a 
and /4) from our best seismic models. Moreover, from this a 
posteriori consistency check, we favour a the more evolved 
seismic models, i.e., those with lower mass and higher core 
overshoot (see Tables|6] and [8ll . 



5.5 Discussion 

We have shown that adopting ?i = 2 for the 
radial mode cannot be exclud ed as was done by 
iDziembowski &i Pamvatnvkhl (|2008l ). and, in fact, explains 
all observed frequencies according to their photometric and 
spectroscopic mode identification except for fg. The iden- 
tifications of the radial order together with the additional 
£ identifications are summarized in Table [T] A perfect 
agreement in frequency values is not achieved. One reason 
might be that slightly different values for {X, Z) than those 
adopted could lead to a better agreement. Another one may 
be that small frequenc y shifts due to rotation may occur a s 
already pointed out by Dziembowski fc Pamvatnvkhl (|2008h , 



An asteroseismic study of the jS Cephei star 12 Lacertae 13 



radial first overtone 



6 p 



3 1 



radial fundamental 



6 p 



4 - 



3 - 



2 - 



1 - 



- 



3 1 



2 
» 



3.5 



4.5 



5.5 6 
Frequency (d"^) 



6.5 



7.5 



Figure 11. Comparison between the theoretical frequencies (m = 0, dots) and the observed frequencies (full and dotted lines) labeled 
with their frequency identifications as listed in Table|5] The full lines correspond to the two axisymmetric modes /2 and /a. The top 
and bottom panels correspond to /2 as the radial first overtone and the radial fundamental, respectively. The theoretical models are 
calculated for X = 0.72, Z = 0.015 and ctov = 0.0. 



although the surface equatorial rotation frequency is only 
a few percent of the measured oscillation frequencies. We 
found evidence for some degree of differential rotation from 
the observed splitting between /i and /a. Using the Ledoux 
splitting, all the other detected frequencies can be well fit- 
ted with one and the same averaged rotation frequency /rot 
inside the star, where /rot is derived from the model fitting 
of /i and /s . In the absence of the observations of multiplets 
we do not have enough constraints to derive the internal ro- 
tation profile as defined in Equation and, consequently, 
to determine its effects on the frequencies and frequency 
splittings. 

iDziembowski fc Pamvatnvkhl l)2008l ) partly included ro- 
tation effects in their modelling and proposed a model with 
the angular rotation velocity more than 4 times higher in the 
centre of the star compared to the surface. They imposed 
that the radial mode is the fundamental without consid- 
ering alternative identifications. Consequently, they could 
not account for an £ — 1 for /y. As our results are the 
same than theirs under the assumption of the radial fun- 
damental mode, we refer to their paper for the properties of 
the resulting seismic model in this case without repeating 
it here (see in particular their Fig. 7 and our Fig. 1111) . Their 



assumption of non-rigid rotation is based on the fact that 
/s cannot be fitted because the closest £ — 2 mode is too far 
to be explained by the splitting due to a uniform rotation. 
We recall that, on the contrary, our models reproduce the 
£ = 1 for /r and the Ledoux splitting for /rot can explain /s 
as the {£,m,n) = (2, 1, —2) mode. As additional argument 
to non-rigid rotation, IDziembowski fc Pamvatnvkhl (|2008l ) 
stated that their model with uniform rotation has no iden- 
tification for fe- We can explain it as the {£ = 4, g2) mode 
or the {£ = 4, gi) mode in agreement with the photometric 
^-identification. Finally, even with a non-uniform rotation, 
they cannot reproduce /s appropriately while they used this 
frequency value to exclude rigid rotation. 

Until now, the only convincing pro of of non-rigid ro- 



tation has been ac hieved for 1^83 6 Cen (|Aerts et al 



Dupret et aLll2004l) and for v Eri (iPamvatnvkh et al.l 



20031: 



2004 



Dziembowski fc Pamvatnvkhl \200i ) . The only other star 
with two observed multiplets is 9 Oph for which the observed 
rotational splittings cannot rule out a rigid rotation model. 
For 12 Lac, we are able to fit the observed frequency spec- 
trum using non-rotating stellar models and Ledoux split- 
tings based on one consistent value /rot for the averaged in- 
ternal rotation frequency derived from the splitting between 
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fi and fs. The deviation between the rotation frequency 
derived from fitting the splitting /i — /a and from the ob- 
served surface rotation frequency points to some degree of 
non-rigidity. With only two components of one and the same 
multiplet observed, we are not able to quantify this state- 
ment. 

Besides our basic stellar modelling, we also checked if 
non-adiabatic computations are able to reproduce the ex- 
citation of the observed frequencies. Indeed, recent works 
revealed shortcomings in the details of the excitation mech- 
anism for t hese massive B-type puls a ting stars. For instance , 
for u Eri (|Pamvatnvkh et al.l [200i : lAusseloos et al.ll2004 ). 
the observed lowest g-mode and highest p-mode frequen- 
cies were not predicted to be excited by standard models. 
The same conclusion is found for 12 Lac. Neither the fre- 
quency fg nor the highest p-mode frequencies (/7 and /g) 
are excited by current models, even with OP opacity ta- 
bles. Since 1991, the agreement between the non-adiabatic 
computations for [3 Cephei models and the observations 
has improved with subsequent update of opacities (see e.g. 
iMiglio et aT]|2007l : lMiglicll2007l '). Consequently, a most prob- 
able explanation for the encountered excitation problem is 
that current opacities are still underestimated in the region 
where the driving of pulsation mod es occurs. We refer to 
iDziembowski fc Pamvatnvkhl (|2008h for an additional dis- 
cussion on this matter. 



6 CONCLUSIONS 

Our study was based on 1820 ground-based, high-resolution, 
high-S/N, multisite spectro scopic measurements spread over 
14 years. As pointed out in [Handler et al.1 (|2006l 'l this effort 
was necessary to identify the previously unknown azimuthal 
order of the pulsation modes and thus to be able to perform 
a detailed seismic modelling of the star. We used the Sim 
4553 A line to derive the pulsation characteristics of 12 Lac. 
In total we find 10 independent frequencies which were also 
detected in photometry. Worth mentioning is that we also 
clearly recover the low-frequency signal. One of our aims 
was to provide a unique identification of as many modes 
as possible for 12 Lac. The important result of combining 
our spectroscopic result s with the ones from the intensive 
photometric campaign (|Handler et al.l l2006l l is the unique 
identification of both land m of the four highest-amplitude 
modes. 

With two state-of-the-art methods, the moment method 
and the FPF method, we were able to identify the azimuthal 
order of the four main modes. We could also give a constraint 
on the azimuthal order of the fifth frequency. Two main 
frequencies were identified as axisymmetric modes. One of 
which has an lvalue of 1 (/s), the other one is a radial 
mode (/2). The two other main frequencies are identified 
as (Zi,mi) — (1, 1) and (Z4,m4) = (2, 1). The conclusion for 
the fifth frequency is that the azimuthal order is likely to be 
positive. In addition, the FPF method could also constrain 
the inclination, i = 48 ± 2°, and the surface equatorial ro- 
tational velocity, Vcq = 49 ± 2kms^^. 

The definite identification of four of the observed modes 
together with some constraints of the wavenumbers (either 
f-values or sign of m) on the other modes allowed us to 
carry out a detailed seismic modelling, under the assump- 



tion of linear oscillation modes. The most significant out- 
come of our modelling is the fact that there is a strong pref- 
erence that the observed radial mode (/2) is the radial first 
overtone. If /2 is taken as the radial fundamental, the fre- 
quency spectrum of 12 Lac cannot be reproduced satisfac- 
torily. We do point out that t he th eory we used, as well as 
IDziembowski fc Pamvatnvkhl (|2008h . does not include non- 
linear effects, while these may be present in 12 Lac (e.g. 
iMathias et"al] 119921 ). It remains to be studied how these 
would effect mode identification and seismic modelling, but 
given the very moderate effect and the limit to only very few 
third-order combination frequencies, we do not expect this 
to alter our conclusions, just as for the case of v Eri, w here 
a linear theory was also used (|Pamvatnvkh et al.ll2004h . 

The complete frequency spectrum of 12 Lac, except the 
low frequency, can be fully identified. Our seismic modelling 
also revealed an excitation problem. Indeed, the range of fre- 
quencies theoretically excited is not large enough compared 
to the observations. This might refiect the fact that opaci- 
ties are still underestimated in the region where the driving 
of pulsation modes occurs. This conclusion is valid for both 
the OPAL and OP opacities. 

Finally, our best-fit models indicate that the overshoot- 
ing parameter Oov has to be lower than 0.4 local pressure 
scale heights to get the best agreement with physical param- 
eters derived from spectroscopic observations. A more re- 
fined seismic modelling requires the computation of a much 
denser grid of models and will be done in future work. 
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